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Abstract 

I present a theoretical discussion of the uncertainties related to the QCD anal- 
ysis of the proton structure function F2{x,Q'^) at small x. The role played by 
the 'unphysical' gluon density is pointed out. It is shown how the study of 
more observables can reduce the theoretical uncertainty and, in particular, an 
alternative method of analysis, based on the introduction of physical anomalous 
dimensions, is suggested. 
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1 Introduction 



One of the main outcome of the physics programme carried out at HERA is the observed 
striking rise of the proton structure function F2{x, Q^) at small values of the Bjorken vari- 
able X (2 ■ 10~^ < X < 10~^) and high values of the momentum transfer (Q^ ^2 GeV^). 

The HERA data on F2 represent the first experimental observation of a cross section 
increasing faster than logarithmically with the energy (see, for instance, Ref. 0). This high- 
energy behaviour in the hard-scattering regime is expected if the underlying dynamics is 
driven by self-interacting massless vector bosons, the gluons. Thus, the steep rise of F2 
certainly confirms one of the basic prediction of perturbative QCD f^. 

However, the main reason why the HERA data have attracted much theoretical at- 
tention goes beyond this point. The issue, indeed, is whether the striking rise of F2 at 
small X calls forth a theoretical interpretation in terms of non-conventional QCD dynam- 
ics. In this context, non-conventional QCD stands for any approach (based either on the 
original BFKL equation ^ or on fc^-factorization [^|13]) in which the small-x behaviour 
of F2{x,Q'^) is studied by resumming logarithmic corrections of the type (as'lnx)" to all 
orders in the strong coupling as- By contrast, no small-x resummation is performed within 
the conventional QCD (or DGLAP @JTl[|) approach: the parton densities of the proton at 
a fixed input scale Ql are evolved in according to the Altarelli-Parisi equation evaluated 
in fixed-order perturbation theory. 

The theoretical motivation for the non-conventional approach based on resummation is 
clear. Since multiple gluon radiation in the final state produces perturbative contributions 
of the type (a^lnx)", as soon as x is sufficiently small (i.e. a^lnl/a; ~ 1), the fixed- 
order expansion in as must become inadequate to describe the QCD dynamics. Thus, in 
principle, the non-conventional approach is certainly more accurate at asymptotically-small 
values of x. The question is whether, in practice, in the HERA kinematic region we are 
already approaching this asymptotic regime. 

In my opinion it is quite difficult to answer this question in the context of the QCD 
analysis of the sole F2. Indeed, the small-x rise of F2 can be obtained as the result of two 
combined effects: the increase of perturbative scaling violation in the small-x region and 
the intrinsic non-perturbative steepness of the gluon density. These perturbative and non- 
perturbative components are mixed up not only on the phenomenological side but, more 
importantly, on theoretical basis. Since the QCD description of a single observable, namely 
F2, requires the introduction of two non-perturbative inputs, quark and gluon densities, the 
distinction between perturbative and non-perturbative components is strongly dependent 
on their own definition rather than on the underlying dynamics. 

A better understanding of QCD physics at small x can be achieved by considering more 
observables and thus (over-) constraining the definition of the parton densities. In partic- 
ular, by simply using two hadronic observables one can formulate the dynamics of scaling 
violation enterely in terms of perturbative quantities that play the role of physical anoma- 
lous dimensions. These anomalous dimensions are unambiguously computable in QCD 
perturbation theory and thus they theoretically appear as golden quantities for comparing 
the conventional and non-conventional approaches. 
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The outline of the paper is as follows. In Sect. § I briefly review the theoretical and 
phenomenological status of the QCD analysis of F2{x,Q'^). In particular, I qualitatively 
discuss the theoretical uncertainties relative to different perturbative approaches and to 
the small-x behaviour of the gluon density. In Sect. ^ I introduce the physical anomalous 
dimensions that control the evolution of F2{x,Q'^) and of the longitudinal structure 
function Fl{x,Q'^). The main features of these anomalous dimensions are discussed in 
Sect. while in Sect. p.2| I present their explicit expressions in resummed perturbation 
theory at small x. Additional observations on the relationship between physical anomalous 
dimensions and parton model are considered in Sect. ^. Section ^ deals with physical 
anomalous dimensions for heavy-flavour structure functions. In particular, it points out 
the kinematical features of the physical anomalous dimensions for observables that depend 
on several large- momentum scales. Some general comments are left to Sect. |^. 



2 The proton structure function F2 and the gluon den- 
sity 



The master equations for the perturbative-QCD study of the proton structure function at 
small X are as follows 



dF^jx.Q^) 
dlnQ^ 



= {e'f)fs{x,Q')+ ... +0{1/Q') , 
= {e}) j\z[Pss{as{Q'),z) fs{x/z,Q' 
+ PsgiasiQ'), z) ~fg (x/z, Q^)]+ ... + 0{l/Q 



(1) 



(2) 



dfg{x,Q'^ 
dlnQ^ 



dz 



Pg,{asm,z)fs{x/z,Q' 



(3) 



where e/ is the electric charge of each quark with flavour /, (e^) = (Z]/ii ej)/A^/ and 
Nf is the number of active flavours. In Eqs. I am using the same notation as in 

Ref. [0. Thus, the singlet density fs and the gluon density fg are related to the usual 
quark (antiquark) and gluon densities fq^ (fq^) and fg by the following relations 



fs{x,Q')=xY^[fq^{x,Q') + fq^{x,Q')] , fg{x,Q')=xfg{x,Q' 
f 



(4) 



and the quark splitting function Pss and Psg are given in terms of the customary Altarelli- 
Parisi splitting functions Pab as follows 

Psg{as, x) = 2NfPq^g{as, x) , Pss{as, x) = J2 [PmA^^s, x) + Pq,qj{as, x)] . (5) 



The dots and the term on the right-hand side of Eqs. respectively denote 

the flavour non-singlet component and higher-twist contributions. The contribution of the 
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non-singlet component to Eqs. (0,0) is quantitatively negligible at small-x and it will be 
neglected throughout the paper. However, its inclusion is formally straightforward (see, for 
instance, Eqs. (5.3-5.5) in Ref. |llQ|). 

The basis for Eqs. ([l|-^ is provided by the factorization theorem of mass singularities 



1^ . According to this theorem the {perturbatively calculable) splitting functions Pabicts, x) 
and the (phenomenological) parton densities fa{x, Q^) are not separately physical observ- 
ables. Only proper combinations (convolutions) of them (for instance, the right-hand 
sides of Eqs. (|1|,3)) are related to measurable quantities. Therefore one has some freedom 
(ambiguity) in defining splitting functions and parton densities. This freedom is called 
factorization-scheme dependence and follows from the fact that hadron scattering cross 
sections cannot be computed within a purely perturbative framework. The factorization 
theorem states that at high momentum transfer Q, the perturbatively non-calculable com- 
ponent of all the cross sections is factorizable in few universal (process independent) parton 
distributions. These parton distributions can be defined using experimental information on 
an equal number of hadronic observables at a certain scale. Having done that, the high-Q^ 
behaviour of all the hadronic cross sections can be unambiguosly (modulo power suppressed 
corrections) computed by using perturbation theory. 

Equations (|I|) and (0) refer to the so-called DIS factorization scheme|^JI3|] . In this 
scheme, Eq. (P actually represents the definition of the singlet-quark density fs- The true 
dynamical information is instead contained in the scaling violations of F2 as described by 
Eq. (^ and by the analogous evolution equation (|^) for the gluon density. 

The Altarelli-Parisi splitting functions entering into Eqs. (0J^) are computable in QCD 
perturbation theory as a power series expansion in as'- 



Pa.(«5,x) = Egj Pir^(a:) , (6) 

and the coefficients Pah~^\'x) in this series can be calculated (at least, in principle) to any 
order n in as- 

Conceptually, the content of the QCD analysis of F2 according to the master equations 
written above is the following. One first computes the splitting functions in Eq. (^ to 
a given perturbative accuracy. Then, by using the experimental information on F2 and 
dF2/d In Q^, from Eqs. (@,^ one can determine the quark and gluon densities as functions of 
X and Q^- Finally, Eq. (H) enters as self-consistency check of the QCD evolution equations. 

In practice, the QCD analysis proceeds as follows. One assigns a certain parametrization 
for the parton densities at a given input scale Qq- Then, inserting this parametrization into 
Eqs. (Ill-H), one can fit the input parameters to the experimental data. 



It is worth emphasizing a point that is independent of the actual procedure used in the 
QCD analysis. Whilst the F2 data uniquely determine the (DIS scheme) quark density, the 
measurement of dF2/d\nQ'^ does not give access directly to the determination of the gluon 

^Equation does not involve (at least, directly) any physical observable. Thus, it takes the same form 
in any scheme. Of course, the gluon splitting functions Pgq and Pgg have to be consistently computed in 
the corresponding factorization scheme. 
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density fg, but rather to that of the product (convolution) Psg ® fg. Since Psg is evaluated 
in a certain theoretical framework (that is, at a given perturbative order or in resummed 
perturbation theory, in a certain factorization scheme and so forth), the ensuing fg turns 
out to be 'theory-dependent'. 



2.1 Fixed-order perturbation theory 

Only the first two terms P^^°'*(a;), Pal\x) of the perturbative expansion (^ of the splitting 
functions are exactly known (i.e. known for any values of x) [1^. In the conventional 
approach these leading order (LO) and next-to-leading order (NLO) terms are used as 
theoretical inputs for the QCD analysis of F2. It turns out that the HERA data can be 
succesfully described []T5| - |T9l by parton densities having the following small-x behaviour 

fs{x,Ql)c^x-''^ , fgix,Ql)c^x-'^ (7) 

with A5 ~ Ag = 0.2 ^ 0.3 at the input scale Qo ~ 4 GeV^. 

Up to the second order in as, the quark splitting functions Pss, Psg in Eq- (0) are 
essentially flat at small x, whilst the gluon splitting functions are steeper and behave as 
follow^ 

P99{(^s,x) ^ -T^ Pgg{as,x) C:^ — , (8) 
Ljp X 

where Ca = Nc, Cp = {N^ — 1)/{2N(), A^'c = 3 is the number of colours and I have defined 
as = CaOs/tt. Thus, the phenomenological success of the NLO QCD approach tells us 
that the rise of F2 at small x is due to the DGLAP evolution in the gluon channel (i.e. it is 
due to Eq. (|^)) combined with a steep behaviour ( ~ ) of the input densities at Qq ~ 
4 GeV2. 



2.2 Resummed perturbation theory 

The basis for the non-conventional QCD approach is provided by the BFKL equation Q]. 
Starting from it, a formalism that is able to combine consistently small-x resummation 
with the QCD factorization theorem has been set up in the last few years. This formalism, 
known as fc^-factorization or high-energy factorization, was first discussed to leading-order 
accuracy in Refs. [^-|^ and then was extended to higher-orders in Refs. In the 

high-energy factorization approach, one ends up with the usual QCD evolution equations 
(namely, Eqs. (|I]-3) in the case of the proton structure function F2) but the splitting func- 
tions Pab{(y.s,x) in Eq. (H) (and, in general, the process-dependent coefficient functions: 
see Sects. are no longer evaluated in fixed-order perturbation theory. They are in- 
deed supplemented with the all-order resummation of the leading (^a§ ln"~^ x), next-to- 
leading (^ag In"""^ x) and, possibly, subdominant (^a^ln^a;, m < n — 2) contributions 
at small x. Note, also, that this resummation can be performed by having full control 



■'■The scheme dependence of the sphtting functions appears only starting from two-loop order. In 
particular, in t'wo-loop order this dependence is pretty mild at small x. 
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of the factorization-scheme dependence of sphtting (and coefficient) functions and parton 
densities [|,|10|,|g. 



The present theoretical status of small-x resummation is the following|. The leading- 
logarithmic (LL) contributions to the gluon splitting functions Pgg{as,x), Pgg{as,x) are 
known 0J^j2^ . Their resummation leads to a very steep (power-like) asymptotic behaviour: 

Pgg{as, x)\asym. ^ Pgg{as , x)\^ym. ~ 0;"^^+^^^ , (9) 

where the power 1 + = 1 + 4Q;5ln2 ~ 1 -|- 2. 65^5 is the so-called intercept of the 
perturbative QCD pomeron. The complete next-to-leading logarithmic (NLL) contribu- 
tions to the gluon splitting functions are not yet known and calculations are in progress 
P^ , |23|p^ . In particular, the contributions proportional to Nf in Pgg have been evaluated 
recently ^4[. Owing to the gluon dominance at high energy, the quark splitting func- 
tions Psg{<ys,x), Pss{c(s,x) do not contain LL contributions. However, the NLL terms are 
completely known |^,[TD| to all orders in perturbation theory. 

Having developed a resummed perturbative expansion to the same (modulo the still 
unknown NLL terms in the gluon sector) degree of theoretical accuracy as the fixed-order 
perturbative expansion, one can set up a fully consistent non-conventional QCD approach 
||T0| . This is accomplished [^,0 by adding leading and next-to-leading logs to one- and 
two-loop contributions (after subtracting the resummed logarithmic terms, in order to 
avoid double counting) in the splitting functions P°'''{as,x), thus obtaining a perturbative 
framework which is everywhere at least as good as the fixed-order expansion, and much 
better becomes small. 



After the first numerical analyses within the /c^-factorization framework, phe- 



nomenological studies based on resummed perturbation theory have been performed during 
the last year [^,^,^. They have shown that, likewise the conventional QCD analysis, the 
non-conventional approach can accomodate the parton densities to provide a description 
of the HERA data on F2. 

The naive explanation for that could be that the inclusion of the resummed logarithmic 
corrections produces a small effect in the Altarelli-Parisi splitting functions. Actually, this 
is not the case. 

Indeed, it is true that the resummation of the leading terms ^agln"~^a; in the gluon 
splitting functions has a moderate impact on the scaling violations of F2 in the kinematical 
range presently investigated at HERA. The situation is however different in the quark chan- 
nel, that is, in the evolution equation (Q). The measured large value of dF2{x,Q'^)/dlnQ'^ 
at small x calls for a quite steep product (convolution) Psg^ fg- In the conventional (fixed- 
order) perturbative analysis this condition can be fulfilled only by choosing a quite steep 
input distribution fg. After resummation of the next-to-leading terms -ag\n^~'^ x, the 
quark sphtting functions Psg{as,x) and Pss{c(s,x) are much steeper than the correspond- 
ing splitting functions evaluated in two-loop orderQ. Thus, the non-conventional approach 

§1 refer to Sects. ^ and |^ for the resummation in the process-dependent coefhcient functions. 



refer to |29| for a more detailed discussion on the smaU-x behaviour of the resummed sphtting 
functions. 
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can succeed in describing the small-x rise of F2 by using parton densities that at the input 



scale Qq are less steep than those needed in the fixed-order approach. 



From this result one may conclude that the conventional and non-conventional ap- 
proaches are phenomenologically equivalent: the HERA data on F2 at sufficiently high 
cannot distinguish steep input densities from steep dynamical evolution. 

The conclusion is instead different from a theoretical viewpoint. Since the resummation 
of the NLL contributions leads to a large effect on F2, there is no justification for truncat- 
ing the QCD perturbative expansion at NLO: the fixed-order expansion approach is thus 
theoretically disfavoured. The only caveat against such a firm conclusion is that the NLL 
terms in the gluon sector are still unknown: they may lead to a large and opposite effect 
with respect to those in the quark channel. 

At the same time, since the known NLL contributions produce large corrections on F2, 
one can expect that subleading terms may still have a sizeable effect. Thus, at present, 
perturbative QCD predictions for the small-x behaviour of F2 suffer from substantial the- 
oretical uncertainties [^. A better understanding of subleading contributions is necessary 



to reduce these uncertainties E5l|30 



The QCD analysis of the sole proton structure function F2, moreover, is affected by 
an even larger indeterminacy related to the difficulty in disentangling perturbative and 
non-perturbative effects. In order to clarify this point, let me brieffy consider the issue of 
the factorization-scheme dependence p0|,p9|,|31|| . 



2.3 Factorization-scheme dependence 

As discussed in the first part of this Section, the parton densities are not physical observ- 
ables and, in particular, they are not calculable in perturbation theory. In the perturbative 
framework they are defined apart from an overall perturbative function. Thus, starting 
from the DIS factorization scheme considered so far, one can introduce a new factorization 
scheme of DIS type| (i.e. a scheme in which the physical identification of the quark density 
with the proton structure function as in Eq. (0) remains valid) by defining a new gluon 
density /j""'") as follows [^9] 



Q') = Q')+ - u{asm. z) f,{x/z, Q') + ^ v{asm. z) fs{x/z, Q' 

Jx Z I Ua 

(10) 

Here u{as,z) and v{as,z) are functions that can be expanded as power series in as and 
vanish for = 0. As for their functional dependence on z, it is quite arbitrary. The 
only constraints are that u{as, z) and v{asi z) contain at most NLL terms of the type 
a5(a5 In z)" for z — > and that these NLL terms are equal in u and v. 

The dynamical evolution equations (§) and (^) can be written in terms of the new 
gluon density ( [T0| ) and of new Altarelli-Parisi splitting functions. The above constraints 
guarantee that the new splitting functions have the same LL behaviour as the DIS-scheme 

"More general factorization schemes, like for instance the MS scheme, are considered in Refs. 31 1. 
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splitting functions, so that the dominant perturbative dynamics is left unchanged. 



The freedom of arbitrarily choosing the factorization scheme is not a particular feature 
of small-x dynamics. The transformation in Eq. (|l^) can be applied in the small-x as well 
as in the large-x regions. In general, its effect amounts to a redefinition of the input parton 
densities that is perturbatively under control. The effect, instead, can be quite large in the 
small-x region because each power of as can be accompanied by an enhancing logarithmic 
factor of In 1/ x. 

In order to quantify the theoretical uncertainty related to the scheme dependence, let 
us consider the simplest case in which the splitting functions in Eqs. (||,||) are evaluated 
to LL accuracy. Thus, we can perform the scheme transformation in Eq. (0) by choosing 
any NLL functions u and v and, in particular, we can set u{as, z) = v{as, z) = Aasz~^°'^ , 
where A and K are constants of order unity. Assuming the extreme case of flat input 
densities, this leads to the following factorization-scheme uncertainty 

= ft'"^{x) - Ux) = I {x-^-^ - 1) ~ ^x-^'^^ . (11) 

This implies that, from the QCD analysis of F2 to LL accuracy, one cannot argue whether 
steep input densities have a non-perturbative origin or rather mimic higher-order pertur- 
bative effects. Of course, using the NLL expressions for the splitting functions one reduces 
the factorization-scheme uncertainty by a factor of as- For the case considered above one 
obtains Sfg{x) ~ asx'^"'^ that, however, still represents a substantial indeterminacy. 

In order to gain more theoretical accuracy one should compute higher perturbative 
orders in the Altarelli-Parisi splitting functions. The same goal can be achieved in a simpler 
manner by eliminating the factorization-scheme uncertainty, that is, by relating the gluon 
density to other physical observables. 

Actually, there is one more reason for studying the small-x behaviour of physical ob- 
servables other than F2. As a matter of fact, the large effect on F2 of the NLL contributions 
in the quark channel might be, in a sense, spurious or, more precisely, related to the use of 
certain factorization schemes pOlEOl. 



For instance (and quite strikingly), one can choose the functions u{as, z) and f (0:5, z) in 
such a way that all the NLL terms in the quark channel are removed from the (new) quark 
splitting functions P^g, Pss and absorbed into the redefinition (^) of the gluon density 
||29|| . The only price one has to pay consists in the introduction of additional NLL terms 
in the gluon splitting functions. It turns out that the effect of these additional terms is 



quantitatively small |26]. Therefore, within this factorization scheme (called SDIS scheme in 
Ref. 1^) the non- conventional approach to NLL accuracy and the conventional approach 
to NLO are, in practice, indistinguishable (apart from the caveat on the unknown NLL 
terms in the gluon channel) as for the QCD analysis of F2. The introduction of the SDIS 
scheme thus provides a more formal argument to explain the phenomenological equivalence 
of the conventional and non-conventional approaches that has been pointed out in Sect. p.2| . 

This equivalence may appear as due to a particular algebraic trick or to fine-tuning 
of the factorization scheme with no physical content. Actually, this is not necessarely 
the case. The factorization-scheme dependence, rather than an ambiguity in higher-order 
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perturbative coefficients, has to be regarded more physically as a parametrization of our 
ignorance in factorizing perturbative from non-perturbative physics. At present, the proton 
structure function F2 (i.e. Eqs. (|ll-^) provides experimental/theoretical information that 
is not sufficiently accurate to disentangle perturbative from non-perturbative dynamics at 
small-x. 

In order to have better control on the perturbative dynamics, one should consider the 
small-x behaviour of other physical observables. Indeed, by means of the transformation 
from the DIS to the SDIS schemes (almost) no trace of small-x perturbative contributions 
is left in F2 and all the resummation effects are moved to other physical quantities. These 
effects can be sizeable. It may also happen that the resummed contributions are almost 
universal, in the sense that, in the kinematic regions that are experimentally accessible, 
they produce very similar quantitative effects in all physical observables. In this case, these 
contributions can be consistently absorbed into the non-perturbative parton densities and 
fixed-order perturbation theory can be safely used throughout. 



3 Factorization-theorem invariants at small x 



The theoretical motivations for studying the small-x behaviour of several different physical 
observables have been pointed out in the previous Section. As discussed in Refs. Il5|j3^ and 
furtherly elaborated on in Ref. this study can be performed by considering properly 



defined i^-factors (ratios of hadronic cross sections), which are factorization-scheme inde- 
pendent. Analogously, one can introduce factorization-scheme invariants that relate the 
scaling violations of different structure functions. These invariants are discussed in the rest 
of this paper. 

Among the observables that one can consider, the longitudinal structure function Fi of 
the proton is becoming increasingly topical. On the experimental side, data on Fl{x, Q"^) at 
small X will be available soon from HERA. On the theoretical side, this quantity is known 
to a sufficient accuracy. 

In order to make more explicit this statement about the theoretical accuracy of Fl, let 
me recall that, using the factorization theorem of mass singularities, F^ is given as follows 

dz 



Fl{x,Q') = (ej)/ - [cl{asiQ'),z) fs{x/z,Q 

•J X Z 



+ Cl{as{Q'),z)~f,(x/z,Q^)\+ ... +0{l/Q^) , (12) 

where fs and fg are the same parton densities that enter into Eqs. (|l]-^ and, as in Eqs. (||,0), 
the dots and the term 0{1/Q'^) respectively denote flavour non-singlet and higher-twist 
contributions. In any given factorization scheme the coefficient functions Cf and Cf in 
Eq. (|l^ are computable in QCD perturbation theory according to the following power 
series expansion 



Cl{as,x) = Y.[^) CT''\x) . (13) 



n=l 



Both the LO and NLO coefficients (x), (x) have been computed for any value 
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of X Correspondingly, in resummed perturbation theory all the NLL terms In" 

in C'l^^\x) are known [|l^]. Owing to this theoretical information, Eqs. (I!]-^ can be sup- 
plemented with Eq. (|T2|) thus eliminating the factorization scheme uncertainty. To this 
purpose one should introduce the parton densities fs and fg extracted from Eqs. into 
Eq. (|l^). Theoretical consistency simply requires that splitting functions and coefficient 
functions are evaluated to the corresponding accuracy, that is, to NLO in the conventional 
approach and including NLL terms in resummed perturbation theory. 



3.1 Physical anomalous dimensions 

The unphysical role played by the parton densities within this context is clear. Indeed, one 
can write down evolution equations that involve only physical observables and perturbative 
quantities. Starting from Eqs. (|I|,|T2|) and performing straightforward algebraic manipula- 
tions, one first express the parton densities fs and fg as functions of F2 and F^. Then, 
inserting the expressions derived in this manner into Eqs. (@J5), one obtains the following 
dynamical equations 



r22(a5(Q'),^) F2 (x/z,Q^) + r2L(as(Q'),^) Fl {x/z,Q' 



d In Jx z 

+ ... +0{l/Q^) , (14) 
dFL{x,Q^) f^dz 



rL2{as{Q^),z) F2 {x/z,Q^) + TLL{as{Q^),z) Fl {x/z,Q' 



d In Jx z 

+ ... +0{1/Q') , (15) 

From a formal viewpoint Eqs. (0,|l5) niay appear equivalent to Eqs. (@,|]). However, 
Eqs. (P^,pr5D relate the scaling violations of two physical observables, namely F2 and Fl, 



to the actual value of the same observables. It follows that the kernels rij{as{Q'^),x) 
(with i,j = 2, L) are physical observables as well. Owing to the formal resemblance to the 
Altarelli-Parisi splitting functions, the kernels Tij{as{Q'^),x) can be considered as physical 
splitting functions. 

The main physical properties of the kernels rij{as{Q'^),x) are that i) each of them is 
consistently computable in QCD perturbation theory (modulo higher-twist corrections that 
are suppressed by some power of l/Q in the hard-scattering regime) and ii) each of them is 
a factorization-theorem invariant, i.e. it does not depend on both the factorization scheme 
and the factorization scale. In other words, from the viewpoint of perturbative QCD, each 
Tij{as{Q'^),x) is completely analogous to the celebrated ratio 

„ a(e^e^ hadrons) , , 

cr(e+e — > ) 

in e~^e~ annihilation. 

The perturbative expansions of the physical splitting functions are the following 



rLL(as,a;) = 2^ — Y\j^'[x^- 



n=l 



27ry ^ ' 27r 



27r 



(17) 
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2 r 



n=A27r 



27r 



rS(x) + gr«w + . 



(18) 
(19) 
(20) 



Note that the expansions for the diagonal kernels Tll, are completely analogous to 
that in Eq. for the Altarelli-Parisi splitting functions. The mismatch in the overall 
power of as between the expansions for the diagonal and non-diagonal {Tl2, ^2l) kernels 
is due to the fact that, from a perturbative viewpoint (or, equivalently, because of the 
validity of the Callan-Gross relation, Fl = 0, in the naive parton model), the longitudinal 
structure function has to be considered as a physical quantity of relative order as with 
respect to F2, i.e. ~ 05^2. Taking this into account, a conventional QCD calculation 
should consistently consider the contributions r^^-* (x) in Eqs. ( p!7| - pO| ) as lowest-order terms, 

r-j'*(x) as next-order terms and so forth. 

Obviously, since the kernels Tij are physical observables, they are renormalization- 
group invariant quantities. It follows that, if computed in fixed-order perturbation theory 
they should exibits the customary dependence on the renormalization scale /i. Thus, to 
be more precise, in the evolution equations (0,0) one has to perform the replacement 



TijiasiQ'^ 



Tiji^asifi"^), Q^Z/i^, x). Equations ( P^pP] ) refer to the perturbative expan- 



sion of Tij for fl = Q . In general one obtains: 



2, Q' 



Tijlasifi ), -^,x 



'as(ji_ 
2n 



2\\P 



r(?fx) 



27r 



(21) 

where 127r/3o = IIC^ ~'^^f is the first coefficient of the QCD /3-function and p = 1 for Tn 
and 122, p = for r2L, p = 2 for rL2- 

This discussion on the perturbative features of the kernels Tij{as-, x) can be summarised 
by saying that they are infrared and collinear safe quantities. Thus, as in the case of the ra- 
tio Re+e- , the x-dependent perturbative coefficients r\j^ are computable by first principles 
starting from parton-level Feynman diagrams and without carrying out any factorization 
procedure of mass singularities. Nonetheless, since higher-order perturbative calculations 
for Altarelli-Parisi splitting functions and process-dependent coefficient functions are al- 
ready available, it is more convenient to relate directly the T^-^^^s to these quantities. 

To the purpose of simplifying the notation it is also useful to introduce the iV-moments. 
For any function g{x), 1 define its A^-moments (7 at in the usual way: 



gj^ = dx x^ ^ g{x) . 
Jo 

Thus, for instance, the evolution equations (|,H) become: 



dF2Mg_ 



(ej) [lss,N{(ys{Q')) fs,N{Q')+lsg,N{asiQ')) Un{Q' 



(22) 



(23) 
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(24) 



where the anomalous dimensions 7afe,Af(tts) are related to the + 1-moments of the 
Aharelh-Parisi sphtting functions, that is, 



lab,N{as)= I dx x^Pab{as,x) = Pab,N+iias) ■ 
Jo 



(25) 



Analogously, the dynamical equations (|1^,|T^) can be rewritten as follows 



rflnQ2 

dFL,Nm 



J- 22, N ^2L,N{as{Q )) Fl,n{Q ) 



rL2,iv(«5(Q'))i^2,iv(g') + J- LL.N («5(Q')) Fl,n{Q' 



(26) 



(27) 



dlnQ2 

where Tij^Nlas) are the physical anomalous dimensions, i.e. the A^-moments of the physical 
splitting functions Tij{as,x). 

The physical anomalous dimensions are related to ■jab,N and to the longitudinal coeffi- 
cient functions in Eq. ([T2|) by the following equations 



LL,N 



c? 



L,N 



1gg,N + lSg,N + 



L,N 



dlnQ^ 



{2i 



J DIS 



L2,N 



CL,N7gq,N - Cf^j^^gg^N + C'f,Af lSS,N - 



lSg,N 



2L,N 



dlnClj^' 
dlnQ2 



7Sg,N 



J DIS 



22, N 



ISS, N — 



^L,N 



r^g -'ySg,N 



-I DIS 



(29) 



(30) 



(31) 



where I have used the shorthand notation r^^Ar = Fj^^ ^^(^^(Q^)), '~iab,N = lab,N{ois{Q'^)), 
^i,N = C'°^(q;s((5^)). In Eqs. ( p^ - |3l| ) the subscript DIS on the right-hand side means 
that the quantities inside the square brackets have to be evaluated in the DIS factorization 
scheme. Obviously, this does not mean that Tj^ jv are scheme dependent. The only point is 
that their expressions in terms of 7afe, n and C£ ^ are more cumbersome if 7afe, at and C£ ^ 
are given in a different factorization scheme. 

Both the Altarelli-Parisi splitting functions and the longitudinal coefficient functions 



are known up to two- loop order. Therefore, using Eqs. (p8| - plD , one can obtain the two 



lowest-order terms vf}, r,-^'* of the physical anomalous dimensions. 
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3.2 Behaviour at small x 



Let me now consider the small-x behaviour of the physical anomalous dimensions. From 
power-counting arguments, it follows that the most singular terms in the perturbative co- 
efficients r\j\x) behave as r\^\x) ~ xP*^"^(x) ~ (Ina;)" or, equivalently, ^[J^n ~ 
in A^-moment spaceQ. As in the case of the Altarelli-Parisi splitting functions, one expects 
two entries in the matrix of the physical splitting functions which contain leading loga- 
rithms. These two entries are those more directly related to the gluon channel and, hence, 
they appear in the evolution equation (^) for the longitudinal structure function. Thus, 
we have r^'^(x) ~ r^'^^(x) ~ (Inx)". The evolution equation ([l^ ) is instead more related to 
the quark dynamics and thus the corresponding anomalous dimensions contain only NLL 
terms, i.e. T2l{^) ~ ^22 {^) ~ (lnx)"~^ 

As in the case of the fixed-order perturbative expansions, Eqs. (p8| - pTD can be used to 
obtain resummed logarithmic expressions at small x for the physical anomalous dimensions. 
The resummation programme carried out in Refs. p|,|10| leads to analytic formulae given in 



terms of the LL contributions to the gluon anomalous dimensions, that is, 

l99,Nias) = 7Jv(«5) + O (asias/Nr) . (32) 

Here, 'jN{cis) is the BFKL anomalous dimension Hj2l[] and is obtained by solving the 
implicit equation 

1 = ^xilNias)) , (33) 

where the characteristic function x{l) is expressed in terms of the Euler ■j/'-function as 
follows 

X(7) = 2^(1) -V'(7)- ^(1-7)- (34) 

Having recalled these results, I am now in a position of presenting all-order resummed 
formulae for the physical anomalous dimensions, starting from the leading components T^l 
and Tl2- 

Since quark splitting functions and longitudinal coefficient functions are subleading at 
small X, Eqs. ( pSj ) and ( p^ ) immediately gives 

TLL,Nicys)=lNicys) +0{asias/N)'') . (35) 

Note also that, using the known next-to- leading results for 75^, and C£, as soon as the next- 
to-leading contributions to 'jgg will be evaluated, one can provide the full NLL corrections 
to Tll [El- 



Incidentally, Eq. ( pSj ) clearly shows that the BFKL anomalous dimensions 7Ar(as'), being 
related to the small-x behaviour of Tn, is a physical quantity. On the contrary, the gluon 
anomalous dimensions 7gg,Ar are factorization-scheme dependent and, in general, one might 
expect that this scheme-dependence affects also their LL behaviour. 



As can be seen from Eqs. ( p9| , |30| , pTD , to the purpose of evaluating the LL terms in rL2, 
as well as the next-to- leading terms in and r22, it is not sufficient to know the leading 

** Note that logarithmic contributions of the type In"^^ x in x-space correspond to multiple poles (1/A^)" 
in A^-space. 
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contributions in the gluon channel, i.e. in and 7^3. One has to use the full information 
provided by the next-to- leading order resummation performed in Ref. This feature 

emphasizes once more that the standard anomalous dimensions and coefficient functions 
are not physical observables. The power counting of small-x logarithms is different for 
physical observables and leading and next-to-leading logarithms in anomalous dimensions 
and coefficient functions get (slightly) mixed up. Without the NLL calculations in Ref. [|lO|] , 
no LL analysis of physical quantities at small x can be carried out. 

Using the results for 73^, •ysg, C*£ obtained in Ref. [jlO[, Eq- (^) gives the following 
expression for the LL terms in 



rL2,N{as) = g - Cf/°)) 7iv(«5) + O {as{as/Nr)] , (36) 



"5 //^CV ^g(O) _ ^S(O)' 

where C^^^ are the A^- moments of the lowest-order coefficient functions C^^^^^ (x) in Eq. ([131) . 
One can see that, eventually, also the small-x resummation in turns out to be propor- 
tional to the BFKL anomalous dimension. As for the NLL terms in Tl2, part of them can 
be obtained from those known for quark anomalous dimensions and coefficient func- 
tions. The remaining terms require the evaluation of the gluon anomalous dimensions jgg 
and 7gg to NLL order and the computation of the DIS-scheme coefficient functions CI to 
next-to-next-to-leading logarithmic (NNLL) accuracy! This feature is consistent with the 
factorization-scheme dependence of Eqs. (0-0). Indeed, it is straightforward to check that, 
after having fixed the factorization scheme to NLL accuracy in Eqs. (|l],|D, only the NLL 
contributions of Pgg in Eq. (^ are unambiguously defined: by properly choosing the scheme 
transformation in Eq. (|T0|), one still has the freedom of arbitrarily defining the NLL terms 
in the non-diagonal gluon splitting function Pgg. In other words, the sole calculation of the 
still unknown gluon anomalous dimensions to NLL order will not be sufficient to provide a 
consistent theoretical framework for the analysis of physical observables to NLL accuracy 
in resummed perturbation theory. 

The evaluation of the next-to-leading contributions to the physical anomalous dimen- 
sions T2L, ^22 enterely relies on the calculations of the quark anomalous dimensions in 
Ref. 1^ and of the longitudinal coefficient functions in Ref. Using these results and 
Eqs. dgjSl), one obtains^] 



"5 -p as 



1 3 , ,■ 

+ o lN[as) 



l-7Ar(a5) 2 



+ 0{aUas/N)') , (37) 



22, N 




r<9 (0) _ r" 



.5(0) 



+ 2 7iv(a5) 



Cf (0) 



- lN[as) 



(38) 



Equations (|37|j38|) provide resummed analytical formulae for T2L and T22 in terms of the 
BFKL anomalous dimension in Eq. (B^). Note that, if one compares the right-hand sides 



^^Note that, consistently with the logarithmic accuracy of the right-hand sides of Eqs. ( ^ , p8| ), the A^- 
moments of the lowest-order anomalous dimensions and coefficient functions can be replaced with their 



values at iV = 0, namely 7^°^ = 0, 7^ 



„(o) 



S(0) 



g,N=0 
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of these equations with the corresponding expressions for 755, 75^, C£ in Ref. |jTO|, one can 



see that Eqs. (|37|j38| ) are remarkably simpler. These equations have to be considered as 
the main scheme-invariant output of the next-to-leading order resummation in the quark 
channel. 

Having presented the main features of the physical anomalous dimensions Tij both 
in fixed-order and in resummed perturbation theory, let me add some comments on the 



dynamical equations ([l^JT5|) 



The first comment regards the theoretical accuracy at small x. Suppose, for instance, 
that the physical anomalous dimensions Tij^isf{as) are evaluated only to LL order in re- 
summed perturbation theory. This implies the following theoretical indeterminacy STn/Tn 
= 0{as{(y.s/NY) or, equivalently, (9(a|(a;5 Ina;)''') in x space. In order to make a direct 
comparison with the discussion in Sect. p.3| , we can parametrize this uncertainty in terms 
of a singular function of the type Aa^x'^""^ . Owing to the convolution structure in 
Eqs. (0,0) and considering the extreme case of flat structure functions, this leads to the 
following theoretical uncertainty 

f^-^a.x--.' , (39) 

Comparing Eqs. ( pT]) and (|39D, we can see that the replacement of unphysical parton 
densities with physical observables (and the ensuing elimination of the factorization-scheme 
dependence) allows one to gain a factor of as in the nominal theoretical accuracy. Of course. 



this is due to the fact that the dynamical evolution equations (0,0) to LL accuracy contain 
more theoretical information than Eqs. to the same accuracy. As a matter of fact, the 
evaluation of the physical anomalous dimensions to LL order is equivalent to the knowledge 
of leading-order splitting functions and next-to-leading order coefficient functions (see the 
discussion above Eq. (|36D). 



Other comments regard phenomenological aspects. The evolution equation (|T5|) for Fl 
is physically analogous to the evolution equation for the gluon density. This analogy is 
particularly clear at small x, because the physical anomalous dimensions Tll and turn 
out to be proportional to the BFKL anomalous dimension. Thus the effects of small-x 
resummation in Eq. (0) can directly be inferred from those studied in Ref. for the 
gluon density. 

The evolution equation (0) for F2 is physically analogous to the evolution equation for 
the quark density. From the expressions in Eqs. (|37| , ^8D we see that the small-x resumma- 
tion effects increase the amount of scaling violation. Equation (|57D, for instance, can be 
rewritten as follows 

+00 

r2L,7v(a5) = l + 2.5 77v(a5) + E(77v(«5))" • (40) 

n=2 

Thus, besides the resummation accomplished by the BFKL anomalous dimension, there 
are further enhancing effects due to the positive definite (although, not large) coefficients 
in the series (^0]). Phenomenological studies of these purely perturbative (i.e. independent 
of the parton densities) effects appear interesting. 
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In general Eqs. ([T^JT5| ) relate measurable values of observables, F2, Fi and their deriva- 
tives with respect to Q^, to perturbative quantities, the physical anomalous dimensions. 
Thus, in the hard scattering regime, these equations provide absolute predictions of per- 
turbative QCD. In practice, the measurement of dFi/dlnQ^ can be quite difficult. In this 
respect, once Fl is measured at a certain value of Q^, from Eq. (^) one can obtain its value 
at any and then one can use Eq. (14) as a test of perturbative QCD that is free from 
non-perturbative parameters. Phenomenological analyses along these lines are in progress 
31. 



4 Parton picture and the unphysical gluon density 



The mathematical steps that are necessary to go from Eqs. (P,p|,|3|JI^) to Eqs. (|ll,|T5) are 
pretty straightforward and, naively, one would be led to conclude that the former equations 
are in one-to-one correspondence with the latter. This is not the case. In order to clarify 
this point let me first consider a case in which a one-to-one correspondence between physical 
and partonic observables can really be estabhshed. 

Suppose we want to evaluate the high-Q^ behaviour of a hadronic observable Fc other 
than, say, F2 and Fl. Suppose also that it is a fiavour-singlet observable measured in lepton- 
hadron scattering processes. Thus, within the partonic framework, we should consider a 
factorization formula analogous to Eq. (|T^). Writing this formula directly in A^-space, we 
have 

Fc,n{Q') = C'cnMQ')) kN{Q') + CI^nMQ')) Un{Q') . (41) 

Using the parton densities as determined from the scaling violations of (for instance) F2 
and Fl, the perturbative QCD prediction for Fq in Eq. (^I]) amounts to the computation 
of two factorization-scheme dependent quantities: the coefficients functions Cq j^{as) and 



Alternatively, we can use Eqs. (PJT^, or Eqs. ([l|,@), to rewrite Eq. (^Tj) as follows 

Fc,n{Q') = Kc2,N{as{Q')) F2MQ') + KcL,NiasiQ')) F^,iv(Q') , 



Fc,n{Q') = Kc2,N{as{Q^)) F2,n{Q') + 



T2C,N{as{Q^)) 

where, using the same notation as in Eqs. (|28| - |30D , we have: 



d\nF2,NiQ^ 
dlnQ2 



C2,N 



{e}) 



L,N 



CL,N 



-I DIS 



(4) L'^i.AfJDIS 



(42) 
(43) 

(44) 



or: 



is: 



C2,N 



S 1SS,N ^g 



lSg,N 



C,N 



2C,N 



DIS 



lSg,N 
^C-, TV J DIS 



(45) 



Equations (|^) and (^) are equivalent to Eq. (|4l|) . The only difference is that the 
factorization-scheme dependence embodied in the coefficient functions and (and in 
the parton densities) has been explicitly eliminated by replacing the parton densities with 
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physical quantities {F2 and in Eq. (|4^ ) or F2 and its Q^-derivative in Eq. (^3D) and in- 
troducing the /^-factors Kc2, Kcl or Kc2, ^2C- These i^'-factors are factorization-scheme 
independent and have the same perturbative properties of the physical anomalous dimen- 
sions (one of them, T2C1 actually coincides with a physical anomalous dimension). As for 
the study of small-x physics and the comparison between fixed-order and resummed per- 
turbation theory, the i^-factors are certainly preferred [^,^ with respect to the coefficient 
functions in Eq. (^T]). However, the perturbative QCD prediction for Fc in Eq. ( ^2]) or (^) 



still involves the computation of two if-factors that replace the two coefficient functions. 

The counting of 'degrees of freedom' is instead different in the case of scaling violations. 
In order to describe the scaling violations of F2 and F^ according to the partonic formulae 
in Eqs. (p],|2|,|3|,p!2|) , one has to assign two input parton densities Ql), fg{x, Ql) (related 
to the low-Q^ behaviour of F2 and F^ or of F2 and its Q^-slope) and to compute sis^ 
quantities in QCD perturbation theory: four flavour-singlet splitting functions Pab{ois,x) 
and two coefficient functions Cl{as,x), Cl{as,x). On the contrary, the solution of the 
dynamical evolution equations (0,|1^) requires two non-perturbative initial conditions at 
the input scale Ql and the calculation of only four quantities, the physical anomalous 
dimensions rij{as,x), in QCD perturbation theory. It is evident that the parton picture 
in Eqs. (|l],|^J|,|l^) introduces spurious perturbative QCD effects. 

Obviously, there is nothing wrong with the parton model or with the (light-cone) Wilson 
expansion for deep-inelastic lepton-hadron scattering. The spurious effects noticed above 
simply follow from the ambiguity in the definition of singlet-quark and gluon densities. 
From a field theory viewpoint, the ambiguity is related to the mixing under renormalization 
of singlet-quark and gluon operators. Owing to the mixing matrix, the renormalization 
prescription has to be specified by four (two, in DIS-type schemes) arbitrary perturbative 
functions. In the partonic framework this ambiguity is unavoidable and ultimately related 
to the fact that no physical current with point-like coupling to gluons does exist. 

The unphysical perturbative contributions that are introduced through the definition 



of the gluon density are responsible for the theoretical uncertainty pointed out in Sect. ^ 
In order to furtherly clarify this aspect, let me discuss another possible effect in resummed 
perturbation theory. 

Suppose that in a certain factorization scheme the resummation of non-leading loga- 
rithmic contributions in the quark or gluon anomalous dimensions produces a singularity 
in the A^-plane at a positive value N = N{as) = CkCtg + . . .. Solving the evolution equation 
(p^) from the input scale Qq to the hard scale Q^, these higher-order terms factorize into 
Qg-dependent and Q^-dependent contributions. In particular, the Qo-dependent factor will 
contain the singularity at = N{as{QQ)) and this singularity will dominate the small-x 
behaviour of F2{x, Q"^) at any unless it is cancelled by a zero in the input parton densi- 
ties fa,N{Qo)- In this case, however, since the parton densities are assumed to be positive 
definite, they must have a singularity at a value of N larger than N{as{Ql)). ^ result 
the small-x behaviour of F2{x, Q^) turns out to be controlled by an unphysical singularity 
that is simply due to the choice of the factorization scheme. 

*This number becomes eight in factorization schemes (like the MS scheme) that are different from the 
DIS scheme. 
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This is certainly an extreme effect. Nonetheless, it shows that the spurious perturbative 
functions that are introduced in the partonic picture may lead to obstructions that cannot 
any longer be removed within the same framework (i.e. without releasing the positivity 
constraint on the parton densities). 



5 Heavy-flavour structure functions 

Most of the discussion in Sect. ^ on physical anomalous dimensions can be repeated for 
other structure functions in deep-inelastic lepton-proton scattering, for instance, the heavy- 
flavour structure functions or F^^. These structure functions are completely analo- 
gous to the customary structure functions F2, Fl with the only additional constraint that 
heavy quarks of mass M are produced in the final state. 

To be precise, all the theoretical formulae in the previous Sections refer to the kinemat- 
ical region ^ and thus neglect corrections of relative order M'^/Q'^. In order to take 
into account the mass effects, one should perform the replacement Fi — > Fi+F^^ . However, 
in the following F^ still denotes the massless contribution to the structure function. 

There are advantages and disadvantages in substituting the heavy-quark structure func- 
tions for Fl in the study of scaling violations. On the experimetal side p5|], the charm con- 



tribution Fg'^ to the small-x behaviour of the proton structure function is certainly more 
easily measurable than Fi. This feature has to be contrasted with a larger theoretical 
uncertainty [^] due to the unknown precise value of the charm mass and with related 
complications considered below. 



In the heavy-flavour case, the analogue of the collinear-factorization formula (jT^) is 
(z = 2,L): 

F?^{i,Q'-,M') = [' ^[c?^'^{as{Q'),^/z;QyM') f2z,Q') 

J ^ z 

+ Cf^''{as{Q^),i/z-Q^/M^)fs{z,Q^)\ . (46) 

Note that F^'^ and the coefficient functions Cf^' " depend on the mass M. Also note that 
in Eq. (^61) I have defined F^^'^ as function of Q^, and the inelasticity variable ^, which 
is related to the customary Bjorken variable x by ^ = x(l + 4M^/Q^). From a theoretical 
viewpoint the scaling variable ^ is preferred to x because it fulfils the kinematical constraint 
< ^ < 1. Thus, considering A^-moments with respect to ^, Eq. (|46|) is diagonalized as 
follows 

F^^{Q'- M') = C^^'^iasiQ'y, Q'/M') f,,^') + C^^^' (asiQ'); Q'/M') fs,N{Q') • 

(47) 

The functions cf^'"" have a perturbative expansion similar to Eq. (p^). There are two 
main differences with respect to the case of C£. For fixed as the heavy-quark coefficient 
functions are not scale invariant because of their explicit dependence on Q"^ /M"^. The ex- 
pansion for Cf^'^as] Q'^/M'^) starts in 0(a|) (i.e. n > 2) and thus only gluons contribute 
to Eq. (^Bf ) at LO (the sensitivity to the gluon density is somehow enhanced). 
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The coefficient functions cf^'^'"'^"' have been fully computed up to NLO (n = 0,1) 
in Ref. ||3^. The corresponding resummed formulae to NLL accuracy^! were obtained in 
Ref. i. 

The evolution equations that involve the physical anomalous dimensions r*^*^ for the 
pair of observables {Fa, F^^^} are the following; 



dF2ix,Q^) 
dlnQ^ 

or, equivalently, in A^-space: 

dF2,Nm 



^dz_ 

X Z 



Tg'^{as{Q'),x/z;QyM')F2iz,Q' 



+ T'^l'{asm,^/^;QVM-') F^'^iz.Q'-M') 

I rl-y r 
5 Z L 

+ TfP{as{Q'),i/z-Q'/M^) F^^{z,Q';M') 



rg%{asm;QVM'')F2,N{Q' 



dings 



= T%%{as{Q'W/M')F,^^{Q^) 
+ T%{asm-,QVM')F^^^{Q'-M') 



dF^9Q(Q^;M' 
d\nQ^ 



(4J 



(49) 



(50) 



(51) 



The relation between the physical anomalous dimensions F'^'^ and the customary split- 
ting and coefficient functions is similar to that in Eqs. (|28| - |3TD for the case of {F2, Fl}, apart 
from the replacement Cl j^ — > Cf '^'"'/ {ej-). Using the same notation as in Eqs. ( ^^ - |3ll) , we 
have: 



r<QQ,s 
r<QQ,g 



rflnCg^ 
rflnQ2 



DIS 



(52) 



i2,N 



{e}) 



/^QQ^ S 



N Igq, N - <^LN Igg, N + (^i^N ISS, N - ^nn a ^^9, N 



r<ww,g 



I rQQ,S ( d In C^^' ^ 



cilnQ2 



DIS 



(53) 



2i,N 



22, N 



lSg,N 

nQQ,g 



DIS 



lSS,N — 



i.N 



DIS 



^The explicit formula for C^'^j^'^ was not reported in Ref. Q and can be found in 



(54) 
(55) 
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5.1 Pert ur bat ive features 



The dynamical evolution equations ( ^ , ^91) (or (|50D,(|5TD) are analogous to Eqs. (|T^|T5|) (or 
(^),(^)). Note, however, a main and important difference: now the physical kernels F'^*^ 
depend not only on as but also on Q"^ /M"^. This dependence is due to the mass- dependence 
of the heavy-flavour coefficients functions C^^'"". In particular, using the identity: 



dlnQ^ 



+ 



(91nQ2 



we can see that the perturbative expansions of the physical anomalous dimensions 

i2 



T?')^ are the following 



rS(«5;QVM^ 



dlnC 



QQ,g{0) 



dlnQ^ 



E 

n=l 



as 

27T 



■.QQ(n-l) 
ii. N 



d\nC- 



QQ,g{0) 



LN 



dlnQ^ 



+ 



as_ 

271 



2tt 



(56) 
' and 



.(57) 



E 

n=l 



2n) 



27r 

2 



F 



QQ(n-l) 
i2,N 



'(0), 



.2,^'(Q7M^) + gFg««(gVM^) + 



Comparing Eqs. ( [57| ) and ([T7|), we see that F^'^ contains a lowest-order contribution (the 



• (58) 



first term on the right-hand side of Eq. (|57|)) that is absent in Tll. This contribution 
leads to kinematical scaling violations, that is, to scaling violations that are independent of 
the running of asiQ"^) and are simply due to the production kinematics of the heavy-quark 
pair. Thus, as for the 'true' dynamical scaling violations, the perturbative QCD calculation 
should consistently consider the coefficients F'^'^*-^'' in Eqs. ( [57|j58| ) as LO terms, F*^*^^^-* as 
NLO terms and so forth. 

Apart from the exphcit (Q2/iVf2)-dependence of the coefficients T'^'^ '^"^~^\Q'^ / M'^) , the 
perturbative expansions of Fg^ and r2'P are completely analogous to those in Eqs. (0) 
and (|20|), respectively. 



Owing to its kinematical origin, the first term on the right-hand side of Eq. (157| ) can 
be eliminated from the physical anomalous dimensions. To this purpouse it is sufficient to 
rescale F^^^ in Eqs. (50,^) by the (factorization-scheme independent) coefficient (j^^'^^^^ 



and, thus, to consider F2^n and Ff^j^/C^ 



QQ,g{0) 



N 



as dynamical variables. Nonetheless, this 



rescaling is not sufficient to exactly put the physical anomalous dimensions F'^'^ on equal 
terms with those in Eqs. ([17| - |20|) . The interplay between kinematical and dynamical scaling 
violations in the heavy-flavour case cannot be avoided beyond the LO. The NLO coefficients 
'^'^^^'i will always depend on the next-to-next-to-leading order (NNLO) (!) 



ii,N ana i at 
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coefficient functions C^^'"'^'^\ Since these have not yet been computed, a fully consistent 



NLO study of the dynamical evolution equation (^91) is, strictly speaking, not feasible at 
present. 

This discussion of the perturbative features of the physical anomalous dimensions 
is not peculiar to the heavy-quark case. It applies to the physical anomalous dimensions 
of any structure function that depends on some other large-momentum scale besides Q^. 



5.2 Small-x resummation 



The power counting of the logarithmic behaviour of F*^*^ at small x is similar to that of the 
physical anomalous dimensions relating F2 and F^. 

The two entries T^'^ and F^*^ have LL contributions. These can be obtained by using 
Eqs. (^,^^ and the known resummed formulae for the heavy-quark coefficient functions 
P] and the quark anomalous dimensions H]. I find: 



rw.^ n2/M2. ^ , , glnifW(7^(a5);gVM^) , , 
Tji, Arias; Q /M ) = 77v(as) H ^nZT^ ^ 0[as{as/N) 



91nQ2 



(59) 



r%%{as;QVM' 



e}) 271 \Ca 



cilnC?g'^(°)(QVM2) 



rflnQ2 



0(as{as/N)'')} , (60) 



where the functions H^^\'y;Q'^/M'^) in Eq. (^) are simply proportional to the i^-factors 
K^^{Q'^ /M"^) introduced in the second paper of Ref. Their eplicit expressions are: 

12 \ 1-7 



4M2 ) 



ih^] 2(1 + 7)7^ + 



■ il + S7,\ F(l-7,l/2;3/2; 



2 + 37 -37^-2(1 + 7) — 



Q2 + 4]^4-2 . 



(61) 



ij(^)(7;QV^') 
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6M2 



7-1 



1 + 



4M2 ) 



F(l-7,l/2;3/2; 



Q2 + 4M2' 



where -F(a, 6; c; z) is the hypergeometric function. 



,(62) 



Owing to the dependence on Q"^ /M"^, the LL behaviour of F^*^ (unlike that of Tn in 
Eq. (P^D) is not simply given by the BFKL anomalous dimension 7Ar(a5). The resummation 
of the LL terms in is achieved through the (Q;s/A^)-dependence of ^N^Ois) and the 7- 
dependence of the function H^'^\'y] Q"^ / Ivp) on the right-hand side of Eq. (p9|) . The LL 
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contributions to F^*^ in Eq. (^) are proportional to F^*^ after subtraction of its lowest- 
order kinematic contribution (cf. Eq. ([57|)). 

The evaluation of F^*^ and F§*^ to NLL accuracy would require the calculation of the 
gluon anomalous dimensions to NLL order and that of the heavy-flavour coefficient 
functions to NNLL order. The NNLL accuracy in C^^"" is demanded by the interplay 
between kinematical and dynamical scaling violations, as discussed in Sect. |5.1|. 



The anomalous dimensions F^*^ and contain only NLL terms at small x. These 
are explicitly given by the following expressions 

7^ ^2i N{(^s; Q /M ) = 2 + 3-fN{as) - 37iv(a5) 



F(l +7jv(a5)) I nr^S/ /Anfc\ ff^o\ 



F(l/2 + 7;v(«5)) Hi^KlNias); Q^M^) 



rQQ /o2/7i^2n '^s j Cf {Q /M ) qq ^2/7i^2n 
^22,n{^s;Q /M ) = — <7^ T^T T2i^^{as;Q /M } 

+ ( i'sIn - % 7i%) I + O [a\{aslNf\ , (64) 



C 



A 



where cq is the heavy-quark electric charge and F(2;) is the Euler F-function. The resumma- 
tion of the logarithmic contributions in Eq. (^) is embodied in the (Q;s'/A^)-dependence of 
BFKL anomalous dimension '^n{(^s) and the 7Ar-dependence of the functions accord- 
ing to Eqs. (^,^). In Eq. (0) the physical anomalous dimension to NLL accuracy 
is expressed in terms of F§*^ through a relation that is analogous to that between F22 and 
F2L in Eqs. 



6 Summary and discussion 

In this contribution I have discussed how the study of different observables can con- 
tribute to our understanding of the dynamics of high-energy hadronic interactions in the 
hard-scattering regime. The main motivation for considering different observables is that 
from the analysis of a single quantity is difficult to disentangle perturbative from non- 
perturbative QCD physics. Of course, we aim to describe both perturbative and non- 
perturbative physics but keeping separate the two aspects can simplify theoretical and 
phenomenological investigations. 

In Sect. |]the interplay between perturbative and non-perturbative dynamics has been 
pointed out in the context of QCD analyses of the small-x behaviour of the proton structure 
function F2{x,Q'^). The factorization theorem of mass singularities provides a representa- 
tion of F2 in terms of phenomenological parton densities and perturbatively computable 
splitting and coefficient functions. As long as the latter have well-behaved perturbative 
expansions, this representation is highly predictive. In the small-x regime, however, higher 
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perturbative orders are strongly enhanced by logarithmic contributions so that, in prin- 
ciple, resummation procedures are mandatory. Thus a physical issue arises: where is the 
boundary between perturbative and non-perturbative phenomena in the hard scattering 
regime? It is quite difficult to tackle this issue by studying the small-x behaviour of the 
sole F2. Indeed, as discussed in Sects. and |2]^, the small-x rise of F2 produced by 
resumming LL and NLL contributions in the perturbative expansion of the splitting func- 
tions is, in many respects (and with the present theoretical and experimental accuracy), 
indistinguishable from a similar rise due to steep parton densities whose evolution is 
performed according to NLO perturbation theory. This uncertainty is formally taken into 
account by the factorization-scheme dependence, as discussed in Sect. ^73| . Owing to this 
dependence, the gluon density may play the role of a hidden variable that, in the case 
of F2, relates different perturbative QCD approaches, namely resummed and fixed-order 
perturbation theory. 

A better theoretical control on perturbative physics can be achieved by exploiting the 
very physical content of the factorization theorem, that is, the universality (process in- 
dependence) of the parton densities. This means that the same parton densities and the 
same perturbative approach have to be used to study the small-x behaviour of different 
physical observables. Universality is particularly evident in the framework of the physical 
anomalous dimensions introduced in Sect. ^. Here I have discussed in detail the case of F2 
and Fi but the method is completely general. 

For any given set fa of parton densities one should consider a set of an equal number 
of hadronic observables Fa- Thus, one can work out the factorization procedure in matrix 
form as follows 

F = Cf, (65) 

where C = Cab is the coefficient function matrix and the simple product structure on the 
right-hand side is usually valid in iV-moment space. Then, it is straightforward to derive 
the following evolution equations 

dF 

TF , (66) 



rflnQ 



where the matrix F of physical anomalous dimensions for the given set of hadronic observ- 
ables is related to C and to the customary matrix 7 of anomalous dimensions as follows 

+ C 7 . (67) 



dlnQ 



While C and 7 are separately factorization-scheme dependent, the physical anomalous 
dimensions (^) are factorization-scheme invariant. As any other infrared and coUinear safe 
observable, they are perturbatively computable apart from corrections that are suppressed 
by some inverse power of Q in the hard-scattering regime. 



In Sect. |3.1| I have considered the physical anomalous dimensions relating the singlet 
components of F2 and F^. These are the most important contributions at small x but the 
physical anomalous dimensions matrix can be introduced in any kinematic region of x. It 
is just sufficient to start from Eq. ( p5| ) by including flavour non-singlet parton densities and 
hadronic observables. 
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In the small-x region the theoretical and phenomenological importance of the evolution 
equations (pBf ) follows from the fact that the small-x perturbative dynamics is completely 
controlled by the physical anomalous dimensions. No spurious perturbative effect (see 
the discussion in Sect. ^ and no subtle interplay between perturbative logarithms and 
steepness of parton densities takes place in Eq. (|66|) . The physical anomalous dimensions 
can be evaluated both in fixed-order perturbation theory and in resummed perturbation 
theory. For any given set of observables and kinematic region of x, one can thus compare the 
two approaches and study the theoretical accuracy of the perturbative expansion. Having 
done that, one can go back to the partonic picture of Eq. (^) and investigate more safely 
the small-x behaviour of the non-perturbative parton densities. 

In Sect. |3.2| , I have presented resummed expressions for the physical anomalous dimen- 
sions Tll, Tl2, ^2L, ^22- Perturbative calculations to NLL accuracy are available for other 
hadronic observables and, in particular, for the heavy-flavour structure functions F^^ and 
F^'^ HI- III Sect. ^ I have considered the corresponding physical anomalous dimensions 
pQQ Using the theoretical approach of Refs. [^0] one can investigate the small-x behaviour 
of other physical anomalous dimensions. In my opinion, some phenomenological studies 
within the framework of physical anomalous dimensions are warranted. 



Acknowledgments. The results reported in this paper were presented and discussed at 
the UK Phenomenology Workshop on HERA Physics, St. John's College, Durham, UK, 
September 1995. I would like to thank Robin Devenish and Mike Whalley for the excellent 
organization of the meeting. 

I also would like to thank Marcello Ciafaloni and Francesco Hautmann for useful comments, 
discussions and a collaboration over many years on these topics. 

Note added. The use of physical anomalous dimensions has been previously advocated 
in the literature (see, for instance, Ref. ||39|) although in different context. I wish to thank 
Georges Grunberg for having pointed out those references to me. 

References 

1. ZEUS Coll., M. Derrick et al, Phys. Lett. 316B (1993) 412, Z. Phys. C65 (1995) 379, 
Z. Phys. C69 (1996) 607, preprint DESY 96-76; HI Coll., I. Abt et al., Nucl. Phys. 
B407 (1993) 515; HI Coll., T. Ahmed et al., Nucl. Phys. B439 (1995) 471. 

2. S. Catani, in Proceedings of Les Rencontres de Physique de La Vallee d'Aoste, ed. 
M. Greco (Editions Frontieres, Gif-sur-Yvette, 1994), p. 227. 

3. A. De Riijula, S.L. Glashow, H.D. Politzer, S.B. Treiman, F. Wilczek and A. Zee, 
Phys. Rev. D 10 (1974) 1649. 

4. L.N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338; E.A. Kuraev, L.N. Lipatov and 
V.S. Fadin, Sov. Phys. JETP 45 (1977) 199; Ya. Balitskii and L.N. Lipatov, Sov. J. 
Nucl. Phys. 28 (1978) 822. 



23 



5. S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B242 (1990) 97, Nucl. Phys. 
B366 (1991) 135, and in Proceedings of the HERA Workshop, eds. W. Buchmiiller 
and G. Ingelman (DESY, Hamburg, 1991), p. 690. 

6. J.C. Colhns and R.K. Elhs, NucL Phys. B360 (1991) 3. 

7. E.M. Levin, M.G. Ryskin, Yu.M. Shabel'skii and A.G. Shuvaev, Sov. J. NucL Phys. 
53 (1991) 657. 

8. S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. 307B (1993) 147. 

9. S. Catani and F. Hautmann, Phys. Lett. 315B (1993) 157. 

10. S. Catani and F. Hautmann, NucL Phys. B427 (1994) 475. 

11. V.N. Gribov and L.N. Lipatov, Sov. J. NucL Phys. 15 (1972) 438, 675; G. Altarelh 
and G. Parisi, NucL Phys. B126 (1977) 298; Yu.L. Dokshitzer, Sov. Phys. JETP 46 
(1977) 641. 

12. See, for instance, J.C. Colhns, D.E. Soper and G. Sterman, in Perturbative Quantum 
Chromodynamics, ed. A.H. Mueller (World Scientific, Singapore, 1989), p. 1 and 
references therein. 

13. G. AhareUi, R.K. EHis and G. Martinelh, Nucl. Phys. B157 (1979) 461. 

14. E.G. Floratos, D.A. Ross and C.T. Sachrajda, Nucl. Phys. B152 (1979) 493; A. 
Gonzalez- Arroyo, C. Lopez and F.J. Yndurain, Nucl. Phys. B153 (1979) 161; G. 
Curci, W. Furmanski and R. Petronzio, Nucl. Phys. B175 (1980) 27; W. Furmanski 
and R. Petronzio, Phys. Lett. 97B (1980) 437; A. Gonzalez-Arroyo and C. Lopez, 
Nucl. Phys. B166 (1980) 429; E.G. Floratos, C. Kounnas and R. Lacaze, Phys. Lett. 
98B (1981) 285. 

15. ZEUS Coh., M. Derrick et al., Phys. Lett. 345B (1995) 576; HI Coll., S. Aid et al., 
Phys. Lett. 354B (1995) 494, preprint DESY 96-39. 

16. M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 

17. CTEQ CoU., H.L. Lai et al., Phys. Rev. D 51 (1995) 4763, preprint MSUHEP-60426 
( |hep-ph/9606399|) . 

18. A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev. D 50 (1994) 6734, Phys. 
Lett. 354B (1995) 155, preprint DTP/96/44 (|hep-ph/9606345|) . 

19. R.D. Ball and S. Forte, Phys. Lett. 335B (1994) 77, Phys. Lett. 336B (1994) 77, Nucl. 
Phys. B (Proc. Supp.) 39BC (1995) 25. 

20. M. Ciafaloni, Phys. Lett. 356B (1995) 74. 

21. T. Jaroszewicz, Phys. Lett. B116 (1982) 291. 



24 



22. L.N. Lipatov and V.S. Fadin, Sov. J. Nucl. Phys. 50 (1989) 712; V.S. Fadin and R. 
Fiore, Phys. Lett. 294B (1992) 286; V.S. Fadin and L.N. Lipatov, Nucl. Phys. B406 
(1993) 259; V.S. Fadin, R. Fiore and A. Quartarolo, Phys. Rev. D 50 (1994) 5893, 
Phys. Rev. D 53 (1996) 2729; V.S. Fadin, Phys. Atom. Nucl. 58 (1995) 1762; V.S. 
Fadin and L.N. Lipatov, preprint DESY-96-020 ( |liep-ph/9602287D; V .S. Fadin, R. 
Fiore and M.I. Kotskii, preprint BUDKERINP-96-35 (|hep-ph/96053571 ). 

23. V. Del Duca, Phys. Rev. D 54 (1996) 989, preprint EDINBURGH 96-3 ( |hep-ph/9604250| ). 

24. G. Camici and M. Ciafaloni, preprint DFF 250-6-96 ( |hep-ph / 9606427]) and to appear. 

25. R.K. Elhs, F. Hautmann and B.R. Webber, Phys. Lett. 348B (1995) 582; F. Haut- 
mann, in QCD and High Energy Hadronic Interactions, Proceedings 30th Rencontres 
de Moriond, ed. J. Tran Thanh Van (Editions Frontieres, Gif-sur-Yvette, 1995), 
p. 133. 

26. R.D. Ball and S. Forte, Phys. Lett. 351B (1995) 313, Phys. Lett. 358B (1995) 365. 

27. A.J. Askew, J. Kwiecinski, A.D. Martin and P.J. Sutton, Phys. Rev. D 47 (1993) 
3775, Phys. Rev. D 49 (1994) 4402; A. J. Askew, K. Golec Biernat, J. Kwiecinski, 
A.D. Martin and P.J. Sutton, Phys. Lett. 325B (1994) 212. 

28. J.R. Forshaw, R.G. Roberts and R.S. Thorne, Phys. Lett. 356B (1995) 79. 

29. S. Catani, Z. Phys. C70 (1996) 263. 

30. S.A. Larin, P. Nogueira, T. van Ritbergen, J. A.M. Vermaseren, preprint NIKHEF- 
96-010 (|hep-ph/9605317| ); S. Forte and R.D. Ball, preprint EDINBURGH-96-14 ( [h^ 



ph/9607291| ), to appear in Proceedings of International Workshop on Deep Inelastic 
Scattering and Related Phenomena (DIS 96), Rome, Italy, April 1996; J. Bliimlein, 
S. Riemersma and A. Vogt, preprint HEPPH-9607329 ( |liep-ph/ 96075^ ), to appear 
in Proceedings of International Workshop on Deep Inelastic Scattering and Related 
Phenomena (DIS 96), Rome, Italy, April 1996. 

31. R.D. Ball and S. Forte, Phys. Lett. 359B (1995) 362. 

32. G. Camici and M. Ciafaloni, Nucl. Phys. B467 (1996) 25. 

33. D.I. Kazakov, A.V. Kotikov, G. Parente, O.A. Sampayo and J. Sanchez- Guilleen, 
Phys. Rev. Lett. 65 (1990) 1535; J. Sanchez- Guilleen, J.L. Miramontes, M. Mira- 
montes, G. Parente and O.A. Sampayo, Nucl. Phys. B353 (1991) 337; E.B. Zijlstra 
and W.L. van Neerven, Nucl. Phys. B383 (1992) 525. 

34. R.S. Thorne, talk at International Workshop on Deep Inelastic Scattering and Related 
Phenomena (DIS 96), Rome, Italy, April 1996; R.G. Roberts and R.S. Thorne, in 
preparation. 

35. K. Daum, talk at International Workshop on Deep Inelastic Scattering and Related 
Phenomena (DIS 96), Rome, Italy, April 1996. 

36. M. Gliick, A. Reya and M. Stratmann, Nucl. Phys. B422 (1994) 37; A. Vogt, preprint 
DESY 96-012 ( |hep-ph/960135^ . 



25 



37. E. Laenen, S. Riemersma, J. Smith and W.L. van Neerven, Nucl. Phys. B392 (1993) 
162; S. Riemersma, J. Smith and W.L. van Neerven, Phys. Lett. 347B (1995) 143. 



38. S. Catani, preprint DFF 254-7-96 ( [hep-ph /96083 10]) , to appear in Proceedings of 
International Workshop on Deep Inelastic Scattering and Related Phenomena (DIS 
96), Rome, Italy, April 1996. 

39. W.A. Bardin and A.J. Buras, Phys. Lett. 86B (1979) 61; E.G. Floratos, C. Kounnas 
and R. Lacaze, Nucl. Phys. B192 (1981) 417; G. Grunberg, Phys. Rev. D 29 (1984) 
2315. 



26 



